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Abstract 



The mesocrystal showing both electrorheological and magnetorheological 
effects is called electro-magnetorheological (EMR) solids. Prediction of the 
overall magnetic properties of the EMR solids is a challenging task due to 
the coexistence of the uniaxially anisotropic behavior and structural transi- 
tion as well as long-range interaction between the suspended particles. To 
consider the uniaxial anisotropy effect, we present an anisotropic Kirkwood- 
Frohlich equation for calculating the effective permeabilities by adopting an 
explicit characteristic spheroid rather than a characteristic sphere used in the 
derivation of the usual Kirkwood-Frohlich equation. Further, by applying an 
Ewald-Kornfeld formulation we are able to investigate the effective permeabil- 
ity by including the structural transition and long-range interaction explicitly. 
Our theory can reduce to the usual Kirkwood-Frohlich equation and Onsager 
equation naturally. To this end, the numerical simulation shows the validity 
of monitoring the structure of EMR solids by detecting their effective perme- 
abilities. 

PACS: 75.30.Gw, 64.70.Kb, 75.50.Kj 
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I. INTRODUCTION 



In 1998 and 1999 a new mesocrystal was reported [1,2] which combines both electrorhe- 
ological (ER) or magnetorheological (MR) effects. This sort of mesocrystal is also called 
electro- and magneto-rheological (EMR) solids. In fact, ER [3] and MR [4] fluids are gener- 
ally particle suspensions in which the particles have large electric polarizability or magnetic 
permeability. In the application of an external electric or magnetic field, the suspended par- 
ticles can form body-centered tetragonal (bet) mesocrystallities, namely ER or MR solids [5]. 
The EMR solid shows very interesting properties when the applied electric field E (in z axis) 
is perpendicular to the magnetic field H (in x or y axis). In detail, in case of dominate elec- 
tric field or dominate magnetic field, EMR solids have thick columns in the dominate field 
direction. These columns have a bet lattice as the ideal structure, too. Recently, a novel 
structural transition in EMR solids was theoretically [1] and experimentally [2] observed 
from bet to face-centered cubic (fee) lattice in the presence of crossed electric and magnetic 
fields as the ratio between the magnetic field and electric field exceeds a minimum value. 

Understanding the magnetic properties of EMR solids is critical to the design of EMR- 
fluid-based devices. Also, these magnetic properties may provide valuable insight into the 
character of the microstructure responsible for their field-dependent rheology as well as 
models of the EMR effect. Since for EMR solids the uniaxial anisotropy occurs naturally, 
the magnetic properties in longitudinal fields (L) should be different from those in tranverse 
fields (T) . Furthermore, the structural transition can affect the effective magnetic properties, 
and the longe-range interaction between the particles (lattice effect) should be expected to 
play an important role as well. Thus, prediction of the overall magnetic properties of EMR 
solids is indeed a challenging task. 

To calculate the effective permeability of EMR solids, the existing methods for cubic 
arrays of spheres [6] or for a suspension containing a dense array of particles [7] can not 
be used directly. Recently, one developed a theory of homogenization to study the effective 
permeability of MR solids with a periodic microstructure [8]. In this paper, we shall present 
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a statistical-mechanical theory, in order to calculate the effective permeability of the EMR 
solids. 

This paper is organized as follows. In Sec. II, by developing the Kirkwood-Frohlich 
equation and using the Ewald-Kornfeld formulation, we present a statistical-mechanical 
theory for the effective permeability of the EMR solids, and the numerical results are given 
as well. This paper ends with a discussion and conclusion in Sec. III. 



II. FORMALISM AND NUMERICAL RESULTS 



A. Contribution of permanent magnetic moments 



For an EMR solid, its effective permeability Jl e is uniaxially anisotropic due to the appli- 
cation of external fields. In detail, the transverse component /4 T ^ (in x or y axis) differs from 
the longitudinal component /4 L ) (in z axis). In this connection, the effective permeability 
Jl e should possess a tensorial form like 
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Since all the permeable particles of the EMR solid have a permanent magnetic dipole moment 
m, it becomes more difficult to derive the expression for /I e . For this purpose, Kirkwood [9] 
and Frohlich [10] introduced a continuum with permeability ju e0 which arises from induced 
magnetization only. Based on it, we shall derive the effective permeability /I e of the EMR 
solid consisting of permeable particles with a permanent magnetic moment m. Let us start 
by seeing each particle with m to have a new tensorial moment m', 

/2 e0 + 2/x 2 i\ 
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and to be embedded in a new host (introduced continuum) of /I e0 , where /i 2 denotes the 
permeability of the nonmagnetic carrier fluid, and / a unit matrix. The denominators in 
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Eq. (2) and in the following equations should be interpreted as inverse matrices. In this 



interactions with the other point dipoles as the particles had, while the magnetizability of 
the particles can be imagined to be smeared out to form a continuum with permeability 
/2 e0 , which will be derived in Sec. (II B). Next, to include the anisotropic feature, we take a 
characteristic spheroid of volume V, which contains n particles. In doing so, the particles in 
the spheroid will be treated explicitly by taking into account the contribution of the particle 
interaction to the effective permeability. In principle, the approximation in this method can 
be made as small as necessary by taking n sufficiently large. Here we should remark that for 
discussing isotropic cases Kirkwood [9] and Frohlich [10] used a characteristic sphere. As a 
matter of fact, no matter for a sphere or a spheroid, each of them should reflect the physical 
properties of the whole suspension. For instance, the number density inside the sphere 
or spheroid should be identical to that of the whole system under consideration. In this 
regard, for the present EMR solid a characteristic sphere is far from being satisfactory, and 
a characteristic spheroid can be used instead so that the uniaxially anisotropic behavior of 
the suspension may be considered more physically. We shall show that the explicit spheroidal 
shape of choice can be determined exactly, see Eq. (19) below. 

All statistical-mechanical theories of the permeability start from 



where B and H denotes the magnetic induction and Maxwell field in the material outside 
the spheroid, respectively. By definition, we write for the magnetization density p t as 



where (M t ) stands for the average total magnetic moment of the spheroid. Here and below 
(• • •) stands for a statistical mechanical average, e.g., 



model each particle is replaced by a point dipole m' having the same non-electrostatic 



B - /i 2 H = Anp t , 
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In this expression, X stands for the set of position and orientation variables of all particles. 
Here U is the energy related to the dipoles in the spheroid, and it consists of three parts: the 
energy of the dipoles in the external field, the magnetostatic interaction energy of the dipoles, 
and the non-magnetostatic interaction energy between the dipoles which are responsible for 
the short-range correlation between orientations and positions of the dipoles. 

Then, the tensorial effective permeability p e of the whole system can be defined as 

BJ = p e U. (5) 

In view of Eqs. (4) and (5), we take one step forward to rewrite Eq. (3) as 

4tt7 

Q2 e -/i 2 /)H = — (M t >. (6) 

Since (M t ) has the same direction as H, it suffices to calculate the average component of 
M f in the direction of H, thus we have 

a-/i 2 7)H=^(M,.e>. (7) 

Here e is the unit vector in the direction of the field. 

In general, p t and (M t ) contain also terms in higher powers of H. Thus, (l/47r)(/2 e0 — 
p?,I)H is the first term in a series development of pil (induced magnetization) in powers of 
H, and must be set equal to the term linear in H of the series development of (M 4 • e)I in 
a Taylor series. So, we obtain 



« 4tt7 fd(M t -e>\ 
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Owing to (M t ■ e) = V(pi + p Q ) and pj = (l/47r)(/2 e0 — yU 2 /)H, we have 



where p Q stands for the orientational magnetization arising from the permanent magnetic 
moments. Rewriting with the external field Hq instead of the Maxwell field H as the 
independent variable we obtain 
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where (M ■ e) = Vp Q . In this case, the external field acting on the spheroid is 



H / 



AT) 



V 










(T) 

M p+(^L M f) )9 (T) 










« (L) 



\ 



(11) 



where the tensorial depolarization factor g represents the spheroid shape. In fact, the 
degree of field-induced anisotropy of the system is determined by how g deviates from 1/3/ 
(isotropic limit). It is worth noting that g will be determined explicitly [see Eq. (19)], and 
that its components satisfy a sum rule 2g^ + g^ = 1 [11]. 
Starting from 
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eventually we have 
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If we use a tensorial Kirkwood correlation factor (3, then we obtain 
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In view of Eqs. (11), (13) and (14), Eq. (10) can be rewritten as 



3kT 



(12) 



(13) 



(14) 



(15) 



where N denotes the number density of the particles. For an isotropic system, namely 
g = 1/3/, Eq. (15) reduces to the usual Kirkwood-Frohlich equation [10] which works for 
permeable particles with a permanent magnetic moment. If g = 1/3/, (3 = 1 and /I e0 = /, 
Eq. (15) reduces to the Onsager equation [12] which treats non-permeable particles with 
a permanent magnetic moment embedded in vacuum. However, it is worth noting that in 
the derivation of the Onsager equation only one particle is considered in the characteristic 
sphere. That is, there is no more correlations between the particle orientations than can be 
accounted for with the help of the continuum method, thus yielding (3 = I . 
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B. Contribution of induced magnetic moments 



Now we are in a position to derive the induced-magnetization-related permeabibility /2 e0 
by performing an Ewald-Kornfeld formulation [13,14] so that the structural transition and 
long-range interaction can be taken into account explicitly. The ground state of the EMR 
solid is a bet (body-centered tetragonal) lattice, which can be regarded as a tetragonal 
lattice, plus a basis of two particles each of which is fixed with an induced point magnetic 
dipole at its center. One of the two particles is located at a corner and the other one at 
the body center of the tetragonal unit cell. Its lattice constants are denoted by a 3 = qi and 
a i(= 0,2) = d-q' 1 / 2 along the z— and x — (y—) axes, respectively. In this case, the uniaxial 
anisotropic axis is directed along z axis. As q varies, the volume of the unit cell keeps 
unchanged, i.e. V c = £ 3 . Thus, the degree of anisotropy of the tetragonal lattice is measured 
by how q is deviated from unity. In particular, q = 0.87358, 1 and 2 1 / 3 represents the bet, 
bec (body-centered cubic) and fee lattice, respectively. 

When one applies an external magnetic field H along x axis, the induced dipole moment 
P are perpendicular to the uniaxial anisotropic axis. Then, the local field H/ c at the lattice 
point R = can be determined. Let us take the transverse component as an example, and 
resort to the Ewald-Kornfeld formulation [13,14] to calculate the local field H ic such that 

2 _ AttP C 2 — C 2 dPn 3 

Hie = P£ £[-7i(^) + x 2 q 2 l2 (R J )} n(G)^exp(-^-) + — |. (16) 
In this equation, 71 and 72 are two coefficients, given by 

= erfc(r/r) + -^exp(-7?V), 

r 5 y/irr 2 

where erfc(?7r) is the complementary error function, and rj an adjustable parameter making 
the summation converge rapidly. In Eq. (16), R and G denote the lattice vector and the 
reciprocal lattice vector, respectively, 

R = l(q~ l l 2 lx + q~ x l 2 my + qnz), 



G = ^-{q ll2 ux + q 1/2 vy + q~ l wz), 

where l,m,n,u,v,w are integers. In addition, Xj and Rj of Eq. (16) are given by, 

Xj = I — ^— — , Rj — \R — ^— — (ax + ay + cz)\, 

and the structure factor 11(G) = 1 + exp[i(u + v + w)/n]. 
So far, let us define a local field factor a, 

a = ^^ L (17) 

It is worth remarking that a is a function of a single variable, namely degree of anisotropy q. 
Also, there is a sum rule 2a^ + = 3 [15]. As q — 1, a — I just represents the isotropic 
limit. Next, we take one step forward to rewrite the well-known Maxwell-Garnett theory for 
isotropic suspensions as [14-16] 

_ f »i-»2 ?> (lg) 



a/2 e0 + (3/ - 5)^ 2 + 
where \i\ stands for the permeability of the particles. This is a developed Maxwell-Garnett 
theory for uniaxially anisotropic suspensions [14,15]. Then it is not difficult to see that the 
depolarization factor g [Eq. (11)] characterizing the shape of the characteristic spheroid of 
choice is determined by 

9 = \a. (19) 
The substitution of Jl e0 [obtained from Eq. (18)] into Eq. (15) leads to /2 e as a result. 



C. Numerical results 

Let us do some numerical simulations. Figure 1 displays and /i^ as a function of q. 
For this figure, we used the Onsager consideration (i.e., assuming the characteristic spheroid 
contains only one particle, (3 = 1), with a focus on the anisotropic effect. As q — 1, this 
system is in the isotropic limit, yielding /4 T ) = (j,( L \ Thus, in Fig. 1 the two points at q — 1 
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are overlapped. It is found that the structural transition of the EMR solid (measured by the 
variation of q) can cause Jl e to change. To some extent, the numerical simulations show the 
validity of monitoring the structure of EMR solids by detecting their effective permeabilities. 

III. DISCUSSION AND CONCLUSION 

Here some comments are in order. An approximation expression for the Kirkwood cor- 
relation factor (3 can be obtained by taking only nearest-neighbors interactions into ac- 
count [17]. In this case, the characteristic spheroid may be shrunk to contain only the i-th 
particle and all the nearest neighbors. It is apparent that (3^ or (3^ will be different from 
1 when there is correlation between the orientations of neighboring particles. When the 
particles tend to direct themselves with parallel permanent magnetic moments, (3^ or (3^ 
will be larger than 1. When the particles prefer an ordering with anti-parallel moments, (3^ 
or (3^ will be smaller than 1. As the EMR solid is subjected to the external magnetic field, 
all the particles can easily direct themselves with parallel permanent magnetic moments. 
In this connection, (3^ or (3^ should be larger than 1, or could approximately equal to 
1 + N c where N c denotes the number of the closest neighboring particles. Nevertheless, once 
N c > 0, the correlation between the nearest particles is included approximately, and this 
(no figures shown here) does not affect the present numerical result on the anisotropic effect 
as (3 {L) = (3 {T) = 1, i.e., N c = 0. In particular, as N c = 4 (i.e., for a BCT lattice), we 
obtain n { p = 17.40 and /4 L) = 18.69. They both are larger than those of N c = due to the 
correlation, as expected. 

The Bergman-Milton spectral representation (BMSR) [18] is an effective method for cal- 
culating the effective dielectric constant of a two-phase composite, and has been successfully 
applied in electrorheological fluids [19], in order to discuss the frequency-dependent complex 
dielectric constant. Alternatively, the BMSR should be expected to work for EMR solids, 
and a favorable comparison between the BMSR and the Ewald summation technique used 
in this work is expected. 
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To sum up, the aim of the present work is to develop a statistical-mechanical theory in 
order to calculate the effective permeability of a new mesocrystal (EMR solid). This theory 
allows one to study the overall magnetic properties of EMR solids, by taking into account 
the anisotropy and structural transition effects and the long-range interaction between the 
suspended particles. Our theory is expected to be of value in computer simulations of 
magnetic/dielectric properties of EMR fluids. 
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FIGURES 

(T) (L) 

FIG. 1. (Color online). (a) n e and (b) He as a function of q. Dot-dashed lines: 
bet (q = 0.87358), bee (q = 1), and fee (q = 2 1 / 3 ). Parameters: m = 2000, /i 2 = 1, 
m = 5.8 x 10~ n emu, / = 0.2, N = 4.2 x 10 6 cm- 3 , and T = 298 K. Solid lines are a guide 
for the eye. 
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